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Abstract — This chapter presents joint interference suppression 
and power allocation algorithms for DS-CDMA and MIMO net- 
works with multiple hops and amplify-and-forward and decode- 
and-forward (DF) protocols. A scheme for joint allocation of 
power levels across the relays and linear interference suppression 
is proposed. We also consider another strategy for joint interfer- 
ence suppression and relay selection that maximizes the diversity 
available in the system. Simulations show that the proposed cross- 
layer optimization algorithms obtain significant gains in capacity 
and performance over existing schemes. 



I. Introduction 

Multiple-antenna wireless communication systems can ex- 
ploit the spatial diversity in wireless channels, mitigating 
the effects of fading and enhancing their performance and 
capacity. Due to the size and cost of mobile terminals, it is 
considered impractical to equip them with multiple antennas. 
However, spatial diversity gains can be obtained when single- 
antenna terminals establish a distributed antenna array via 
cooperation |[T|- |[3|. The use of cooperative strategies can 
lead to several types of gains ||4|, lITSl . namely, pathloss, 
diversity and multiplexing gains. Pathloss gains allow a sig- 
nificant reduction in the transmitted power for an equivalent 
performance, can increase the coverage llTSi and enhance the 
interference suppression capability H, llT3l . The diversity 
gains improve the performance of the wireless system with 
respect to the probability of error because the transmission 
of multiple copies of the signals reduce the probability that 
the message will not be received correctly. The multiplexing 
gains IIT4I . which correspond to the additional number of bits 
that the system can transmit as compared to a single-antenna 
link, can be obtained when a designer can use relays to form 
independent channels and increase the rate of communication. 

Despite the many advantages in terms of gains as previ- 
ously outlined, cooperative communications also entail some 
disadvantages such as signalling overheads |fT3T|, more com- 
putationally complex scheduling algorithms m and increased 
latency ifTSI . For this reason, it is important to weigh the pros 
and cons of cooperative techniques prior to their adoption and 
consider the practical scenarios of interest IITtI . Motivated by 
their performance and diversity gains, cooperative techniques 
are now being considered for the next generation of mobile 
networks lfT2l . [18|, |19|. In cooperative systems, terminals 
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or users relay signals to each other in order to propagate 
redundant copies of the same signals to the destination user or 
terminal. To this end, the designer must resort to a cooperation 
protocol such as amplify-and-forward (AF) ||3], decode-and- 
forward (DF) |3|, |20| and compress-and-forward (CF) [1211 . 

In order to obtain the benefits of cooperative techniques, 
designers must address a number of problems that are encoun- 
tered in cooperative wireless systems. These problems include 
physical-layer strategies such as synchronization, interference 
mitigation, and parameter estimation. However, designers also 
have to consider a number of associated problems that be- 
long to higher protocol layers and include the allocation of 
resources such as power, relays and rate. These tasks present 
an opportunity to perform cross-layer design and to obtain very 
significant gains in performance and capacity for cooperative 
wireless networks. This chapter is concerned with cross- 
layer design techniques for cooperative wireless networks and 
investigates the benefits of approaches that jointly mitigate 
interference and perform resource allocation. 

In this chapter, we will consider two types of schemes, 
namely, direct-sequence code-division multiple access (DS- 
CDMA) Q, H and multi-input multi-output (MIMO) 0, (6| 
systems. The former is of fundamental importance in wireless 
ad-hoc and sensor networks I?), whereas the latter is one of the 
main ingredients of future wireless cellular networks. When 
implementing cooperative techniques in wireless systems, de- 
signers often consider the transmission technologies available 
and their suitability to certain applications. Therefore, the 
concept of distabuted antenna arrays can be easily extended 
to techniques such as MIMO 0, (§1 and DS-CDMA systems 
Q, 10. 

In the context of MIMO systems, one can obtain substantial 
multiplexing 15], 161, ifTTl and diversity gains ||9l, ifTOl with 
the deployment of multiple antennas at both ends of the 
wireless system. MIMO technology is poised to equip most 
of the future wireless systems and can be incorporated in 
conjunction with other transmission systems. There are two 
basic configurations which exploit the nature of the wireless 
channel: spatial multiplexing IfTTl and diversity IfTOl . Spatial 
multiplexing relies on the concept of forming individual data 
stream between pais of transmit and receive antennas. The 
capacity gains of spatial multiplexing grow linearly with the 
minimum number of transmit and receive antennas ||5], IS 
and allow a MIMO system to obtain a considerable increase 
in data rates. Diversity configurations adopt space-time codes 
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[|9l , IfTOl to transmit data from the antennas at the transmitter 
and can obtain a lower probabiHty of error 

DS-CDMA systems are a key multiple access technology 
for current and future wireless communication systems. Such 
systems rely on the idea of transmitting data with the aid of 
unique signatures, which are also known as spreading codes. 
These signatures are responsible for spreading the informa- 
tion in frequency, and allow the system to have multiple 
users on the same channel. The advantages of DS-CDMA 
include good performance in multi-path channels, flexibility 
in the allocation of channels, increased capacity in bursty 
and fading environments and the ability to share bandwidth 
with narrowband communication systems without deteriora- 
tion of cither's systems performance Q, iS]. Demodulating 
a desired user in a DS-CDMA network requires processing 
the received signal in order to mitigate different types of 
interference, namely, narrowband interference (NBl), multi- 
access interference (MAI), inter-symbol interference (ISI) and 
the noise at the receiver The major source of interference in 
most CDMA systems is MAI, which arises due to the fact 
that users communicate through the same physical channel 
with non-orthogonal signals. 

The similarities between MIMO and CDMA systems in- 
clude their mathematically similar descriptions and their fun- 
damental need for interference mitigation. Indeed, the data 
streams of MIMO systems operating in a spatial multiplexing 
configuration are equivalent to the users of a DS-CDMA 
system. In order to separate data streams or users, a designer 
must resort to detection techniques Il22l . which are very similar 
when applied to either MIMO or DS-CDMA. The optimal 
maximum likelihood (ML) detector is often too complex to be 
implemented for systems with a large number of antennas. For 
this reason, designers often resort to suboptimal solutions that 
an attractive trade-off between performance and complexity. 
These include the sphere decoder (SD) algorithms ||231 . hnear 
detectors 1221 . the successive interference cancellation (SIC) 
approach ifTTl . the parallel interference cancellation (PIC) 
llMl and the decision feedback (DF) detectors |[39l, 1691 
are techniques that can offer an attractive trade-off between 
performance and complexity. These detection algorithms can 
be combined with cross-layer design techniques for enhanced 
interference mitigation and improved overall performance. In 
this chapter, we are specifically interested in exploring the 
advantages of linear detection with power allocation, data 
stream and relay selection. 

A. Prior and Related Work 

Prior work on cross-layer design for cooperative and multi- 
hop communications has considered the problem of resource 
allocation 1241 . l25l in generic networks. These include power 
and rate allocation strategies. Related work on cooperative 
multiuser DS-CDMA networks has focused on the assessment 
of the impact of multiple access interference (MAI) and 
intersymbol interference (ISI), the problem of partner selection 
II20I . l26l . the bit error ratio (BER) and outage performance 
analysis l27]| . and training -based joint power allocation and 
interference mitigation strategies ll28l . ||29l . Previous works 
have also considered the problem of antenna selection, relay 
selection (RS) and diversity maximization, which are central 



themes in the MIMO relaying literature ||3T| - ||33l . However, 
current approaches are often limited to stationary, single relay 
systems and channels which assume the direct path from the 
source to the destination is negligible l32l . 

Most of these resource allocation and interference miti- 
gation strategies require a higher computational cost to im- 
plement the power allocation and a significant amount of 
signalling, decreasing the spectral efficiency of cooperative 
networks. This problem is central to ad-hoc and sensor net- 
works II30I that employ spread spectrum systems and require 
multiple hops to communicate with nodes that are far from 
the source node. This is also of paramount importance in 
cooperative cellular networks. 

B. Contributions 

In this chapter, we present joint interference suppression 
and power allocation algorithms for DS-CDMA and MIMO 
networks with multiple hops and AF and DF protocols. A 
scheme that jointly considers the power allocation across the 
relays subject to group-based power constraints and the design 
of linear receivers for interference suppression is proposed. 
The idea of a group-based power allocation constraint is 
shown to yield close to optimal performance, while keeping 
the signalling and complexity requirements low. A constrained 
minimum mean-squared error (MMSE) design for the receive 
filters and the power allocation vectors is developed along with 
an MMSE channel estimator for the cooperative system under 
consideration. The linear MMSE receiver design is adopted 
due to its mathematical tractability and good performance. 
However, the incorporation of more sophisticated detection 
strategies including interference cancellation with iterative de- 
coding ||39l and advanced parameter estimation methods P9l 
is also possible. In order to solve the proposed optimization 
problem efficiently, a method to form an effective group of 
users and an alternating optimization strategy are presented 
with recursive alternating least squares (RALS) algorithms for 
estimating the parameters of the receiver, the power allocation 
and the channels. A joint relay selection and transmit diversity 
selection strategy for MIMO networks with linear receivers 
is also proposed which optimizes relay transmissions with 
minimal feedback requirements. Effectively a novel approach 
to 1-bit power allocation, two joint discrete optimization 
functions are formed which are solved using discrete stochastic 
algorithms. 

C. Organisation of the Chapter 

The chapter is organized as follows. Section II describes 
cooperative DS-CDMA and MIMO system models with mul- 
tiple hops. Section III formulates the problem, details the 
constrained MMSE design of the receive filters and the power 
allocation vectors subject to a group-based power allocation 
constraint, and describes an MMSE channel estimator An 
extension to cooperative MIMO systems is also presented and 
discrete optimization problems are formulated to jointly select 
the optimal relays and their transmit antennas. Section IV 
presents an algorithm to form the group and the alternating 
optimization strategy along with RLS-type algorithms for 
estimating the parameters of the receiver, the power allocation 
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and the channels. For the solution of the combinatorial prob- 
lems posed by the relay selection strategy, a pair of discrete 
stochastic algorithms are introduced and their joint operation 
detailed. Section V presents and discusses the simulation 
results and Section VI draws the conclusions of this work. 



II. System and Data Models of Cooperative 
Wireless Systems 

In this section, we consider system and data models of 
cooperative wireless systems. The basic idea is to use a 
linear algebra approach to describe models of the cooperative 
systems of interest. In particular, we focus on DS-CDMA 
and MIMO systems and we present a unified approach to the 
description of these systems. 



A. Cooperative DS-CDMA System and Data Model 



respectively. 
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Fig. 1. (a) Uplink and (b) downlink of the cooperative DS-CDMA system. 



Let us first consider a synchronous DS-CDMA network 
with multipath channels. The DS-CDMA system operates with 
QPSK modulation, K users, N chips per symbol and L as 
the maximum number of propagation paths for each link. 
An outline of the system is depicted in ([T]i. The system is 
equipped with AF and DF protocols that allow communication 
in multiple hops using n,- fixed relays in a repetitive fashion. 
We assume that the source node or terminal transmits data 
organized in packets with P symbols, there is enough control 
data to coordinate transmissions and cooperation, and the 
linear receivers at the relay and destination terminals are 
synchronized with their desired signals. The received signals 
are filtered by a matched filter, sampled at chip rate and 
organized into M x 1 vectors r^d, fgrj and r^^d, which 
describe the signal received from the source to the destination, 
the source to the relays, and the relays to the destination. 
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J2°'r,dCkhr^d,kbk+Vr,d 
k=l 

+ rirjd, 

j = l,...,nr, i = l,...,P 

where AI = N + L — 1, P is the number of packet symbols, 
Tir + l is the number of transmission phases or hops, and 
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rir is the number of relays. The vectors Ugd, ngr^ and rir-d 
are zero mean complex Gaussian vectors with variance 
generated at the receivers of the destination and the relays from 
different links, and the vectors rj^^, rf^^ and 77,. represent the 
intersymbol interference (ISI). The amplitudes of the source to 
destination, source to relay and relay to destination links for 
user k are denoted by a^^, a^^^ and a^.^i, respectively. The 
quantities bk and hk represent the original and reconstructed 
symbols by the AF or DF protocol at the relays, respectively. 
The M X L matrix Ck contains versions of the signature 
sequences of each user shifted down by one position at each 
column as described by 
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where Ck = [cfe(l), Cfc(2), 



Ck{N)] stands for the 



signature sequence of user k, the L x 1 channel vectors from 
source to destination, source to relay, and relay to destination 
are hgd.k, hsr^.k, hp.j^ k, respectively. By collecting the data 
vectors in Q (including the links from relays to the destina- 
tion) into a J X 1 received vector at the destination, where 
J = {rir + 1)M, we obtain 
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Rewriting the above signals in a compact form yields 

K 

r[i] = ^ Bk\i\Ak\i\ Ckhk\i\ +r][i\ + n[i] 
PfcW 



k=l 
K 



k=l 
K 

= ^Pk\i\Bk\i\ak\i\+'n[i]+ 71% 

k=l 

where the J x (ji^ + 1)L matrix Ck = diagjCfe . . . C^} 
contains copies of Ck shifted down by M positions for each 
group of L columns and zeros elsewhere. The Q x 1 vector 
hk[i], where Q = (n^ + 1)L contains the channel gains of 
the links between the source, the relays and the destination, 
and Pi.[i] is the effective signature for user k. The {ur + 1) x 
(n^ + l) diagonal matrix Bk[i] = diag(&fc[j] • • • bl"'^[i]) 

contains the symbols transmitted from the source to the 
destination {bk[i]) and the n,. symbols transmitted from 

the 



the relays to the destination . . . on 

main diagonal, and the J x J diagonal matrix -Bfc[?'] 
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[i]^lM), where (g) 
denotes the Kronecker product and Im is an identity matrix 
with dimension Af . The {ur + 1) x 1 power allocation vector 
ak[i] = [a^'^ aj!^^ . . . a^^ j]"^ has the amplitudes of the links, 
the {rir + 1) X (rir + l) diagonal matrix Ak[i] is given by 
Ak[i] = diag{afe[i]}, and the J x J diagonal matrix Ak[i] = 
Wsd <d (g) Ia/ . . . a^„^rf (g) ImV- The J x (n, + 1) 
matrix Pk has copies of the effective signature shifted 
down by M positions for each column and zeros elsewhere. 
The J X 1 vector ri[i] represents the ISI terms and the J x 1 
vector n[i] has the noise components. 



B. Cooperative MIMO System and Data Model 

Let us now consider a synchronous MIMO system model, 
which has similarities with the DS-CDMA system model of 
the previous subsection. We consider a narrowband MIMO 
system with flat fading channels, QPSK modulation, K trans- 
mit antennas, and AI receive antennas as illustrated in Fig. 
121 The cooperative MIMO network is equipped with DF 
protocol that allows communication in = 2 hops using 
rir fixed relays in a repetitive fashion where a non-negligible, 
direct source to destination link exists during the first phase. 
We assume that the source node or terminal transmits data 
organized in packets with P symbols, there is enough control 
data to coordinate transmissions and cooperation, and the 
linear receivers at the relay and destination terminals are 
synchronized with their desired signals. It should be noted 
that the MIMO and CDMA system and data models are 
mathematically equivalent and the main difference is that 
we employ for the MIMO version a spreading code matrix 
Ck = l. 

The received signals are filtered by a matched filter, sampled 
at chip rate and organized into M x 1 vectors Tgd, fsr, and 
r^d, which describe the signal received from the source to 
the destination, the source to the relays, and the relays to the 



destination, respectively, 
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where P is the number of packet symbols, = 2 is the 
number of transmission phases or hops, and rir is the number 
of relays. The vectors Usd, Wsr and rir d are zero mean 
complex Gaussian vectors with variance generated at the 
receivers of the destination and the relays from different links, 
and the vectors rj^^, rj^r and ry^ ^ represent the intersymbol 
interference (ISI). 

The amplitudes of the source to destination, source to relay 
and relay to destination links for user k are denoted by a^"^, 
a'^r and aj!.^, respectively. The quantities bk[i] and bk[i] 
represent the original and reconstructed symbols by the AF 
or DF protocol at the relays, respectively. The Af x 1 spatial 
channel vectors from source to destination, source to relay, 
and relay to destination are hsd,k, hgr^^k, hr-d,k, respectively. 
By collecting the data vectors in (|5]l (including the links from 
relays to the destination) into a J x 1 received vector at the 
destination, where J = UpM for MIMO systems, we obtain 
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Fig. 2. Block diagram of 2-pliase cooperative MIMO system. 
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Rewriting the above signals in a compact form yields 

r[i\ = ^ Bk\i\Ak[i\hk[i] + n\i\ 

= ^Pk\i\Bk\i\ak\i\+n[i], 

k=l 

The J X 1 vector hk[i\ contains the spatial channel gains of 
the Unks between the source, the relays and the destination. 
The rip x rip diagonal matrix Bk[i\ = diag(6fc[i] b]!^'''[i\), 
for 1 < fc < K, and Bk[i\ = diag(0 bl^'^[i\), for 
K < k < Krir, contains the symbols transmitted from 
the source to the destination and the 77,,. symbols 

transmitted from the relays to the destination {h]^'''[i]) on 
the main diagonal, and the J x J diagonal matrix Bk[i] = 
diag(5fe[z]0JM K^'IA^Im), for 1 < fc < K, and 
Bk[i] = diag(0{g)/A/ bl^%]®lM). for K < k < Kur,. 
The rip X 1 power allocation vector ak[i] = [a^'^ . . . arkd]"^ 
has the amplitudes of the links, the rip x rip diagonal matrix 
Ak[i] is given by Ak[i] ~ diag{afc[i]}, and the JxJ diagonal 
matrix Ak\i] = [a'^^^lM ■ • ■ a^.^^^^ I m]'^ ■ The J x Up 
matrix Pk has copies of the spatial signatures hk[i] shifted 
down by M positions for each column and zeros elsewhere. 
The J X 1 vector n[i] has the noise components. 

III. Joint MMSE Receiver Design, Power 
Allocation, Relay Selection and Channel 
Estimation 

In this section, our aim is to describe techniques to mitigate 
interference and allocate the power and select the best relays 
according to the mean square error (MSB) criterion. We 
present a joint receiver design and power allocation strat- 
egy using constrained linear MMSE estimation and group- 
based power constraints along with a linear MMSE channel 
estimator The interesting aspect of the group-based power 
constraints is that a designer can choose a subset of users 
or data streams for power adjustment. Another technique that 
is detailed here is a method called transmit diversity selection 
(TDS) which operates with the linear MMSE receivers. With 
the TDS technique the relay selection is used to jointly 
optimize the selection of antennas in a strategy equivalent to 
a 1-bit transmit antenna power allocation. 

In order to describe the techniques necessary for inter- 
ference mitigation and resource allocation, we introduce an 
alternative way of expressing the J x 1 received vector in (|7]i. 
Our goal is to separate the subset of users or data streams 
that will be used for resource allocation with the group-based 
power constraints. The modified J x 1 received vector can be 
expressed as 

r\i] = Ps [i]Bs [i]as,k M + W-Bfe [i]ak [i] + viA + " W , 

(8) 

where <S = {.Si, 1S2, . . . , <Sg} denotes the group of 
G users to consider in the design. The J x G{nr + 
1) matrix Pg = [Ps^ Ps-2 ■ ■ ■ Psg\ contains 
the G effective signatures of the group of users. The 

G{nr + 1) X G{nr + 1) diagonal matrix Bs\i] = 
diag(65jr] %\'[i]...~bl-,'[^ ... bsa[i] %'^[i]---bZ'm 



contains the symbols transmitted from the sources to 
the destination and from the relays to the destination 
of the G users in the group on the main diagonal, 
the G{nr + 1) x 1 power allocation vector 05^(1] = 

[«S W • ■ • <IM • ■ • > «f." W <m ■ ■ ■ <% of 

the amplitudes of the links used by the G users or data streams 
in the group. 

A. Linear MMSE Receiver Design and Power Allocation 
Scheme with Group-Based Constraints 

The linear MMSE interference mitigation for user or 
data stream fc is performed by the receive filter Wk[i] ~ 
• ■ • , Wfc^j[z]] with J coefficients on the received data 
vector r[i] and yields 

Zk[i]=w^[i]r[i], (9) 

where Zk [i] is an estimate of the symbols, which are processed 
by a sheer Q{-) that performs detection and obtains the desired 
symbol as bk[i] = Q{zk[i])- 

Let us now detail the linear MMSE-based design of the 
receivers for user or data stream k represented by Wk [i] and for 
the computation of the G{nr + 1) x 1 power allocation vector 
o.s.k[i\- This problem can be cast as the following constrained 
optimization 

K"*. «S] = arg min E[i\bk[i] - t^f [i]r[i]n 
subject to af fe[i]a5_fc[i] = Pq, 

(10) 

In order to obtain expressions for the receive filter Wk [i] and 
the power allocation vector a5.fe[i] subject to the group-based 
power constraints, we need the help of the method of Lagrange 
multipliers ( fTOl l ifSTI that transforms a constrained optimization 
into an unconstrained one. The MMSE expressions for Wk [i] 
and a5 jj[r] are given by 

Ck ^ E[\bk[i\ - w^[t\{Ps\i]Bs\i]as,k[i\ 

+ HBfc [i]a, [{] + r,[i] + n[{\) f] j) 

k^S 

+ Afc(as,fc[z] - Pg), 

where Afe is a Lagrange multiplier 

Since the Lagrangian in (fTTT i is a function of both Wk[i] 
and a5.fe[r], we need to employ a strategy for optimization 
the function with respect to both parameter vectors. The main 
idea is to fix one of the parameter vectors and compute the 
gradient terms with respect to the other parameter vector 
that minimizes the Lagrangian and obtain the expression of 
interest. In particular, an expression for 05 is obtained by 
fixing iffc[i], taking the gradient terms of the Lagrangian and 
equating them to zero, which yields 

as.k[i] = (-R5,feW + Xkiy^PsAi] (12) 

where the G(ri,. + 1) x G{nr + 1) covariance matrix 
Rs.ki^ = E[B^[i]Ps[^]wk[i]w^[i]Ps[i\Bs[i]] and the 
vector = E[bk[i]B^[i]P^[i]wk[i]] is a G(n, + 1) x 1 

cross-correlation vector The Lagrange multiplier Afe plays 
the role of a regularization term and has to be determined 
numerically due to the difficulty of evaluating its expression. 

In order to compute the expression for t(;fe[i], we fix ^[i], 
calculate the gradient terms of the Lagrangian and equate them 
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to zero which leads to 



Wk[i\ = R ^WPfcW, 



(13) 



where the covariance matrix of the received vector is given 
by R[i\ = E[r[i\r"[i\] and Pf.[i\ = E[bl[i\r[i\] is a J x 1 
cross-correlation vector The quantities R[i] and depend 
on the power allocation vector a5 fc[i]. The expressions in 
( fT2] i and ( flSb do not have a closed-form solution as they have 
a dependence on each other Moreover, the expressions also 
require the estimation of the channel vector hk[i]- Thus, it 
is necessary to iterate (fTZt and ( fTJt with initial values to 
obtain a solution and to estimate the channel. The network 
has to convey the information from the group of users which 
is necessary to compute the group-based power allocation 
including the filter tOfc[i]. The expressions in ( fT2l i and ( fT3T l 
require matrix inversions with cubic complexity ( 0{{J)^) 
and 0((GK + 1))^)- 

B. Transmit Diversity Selection and Relay Selection 

In this subsection, we explore the idea of transmit diversity 
selection (TDS) and relay selection (RS) and how they can 
be used to improve the performance of cooperative systems. 
In cooperative wireless systems with multiple relays, there 
are links that have very poor propagation conditions that can 
degrade the performance of the overall system. These links can 
be identified and removed from the operation of the system 
via TDS and RS. To this end, we formulate a TDS and RS 
strategy for a 2 DF MIMO network as a discrete combinatorial 
MSB problem which optimizes the use of the channels of 
the second phase via 1-bit power allocation |551. It turns out 
that the problems of TDS and RS are combinatorial problems 
which require either an exhaustive search or some relaxation 
approach. We specify that a subset of Ksub antennas of the 
rirK relay antennas are active at each time instant in order 
to reduce the optimization complexity but also to ensure a 
minimum available level of diversity. The destination node's 
MSB TDS optimization function is given by 

r - (14) 



arg mm 



where T, = diag(ai^^ . . . a,^^, . 



-.K 



i) and 



= {0,1}, and Wk is the linear MMSB filter for the 
l^th sy[jj]3ol. Under the assumption of no inter-relay com- 
munication and that each data stream is allocated to its 
correspondingly numbered transmit antenna at each relay, the 
set fix has a cardinality of {flrl — and contains all 

possible combinations of relay transmit antennas patterns. The 
performance and complexity of solutions to (fl4l i depend on 
{flrl and its elements. However, I^ItI is significant even for 
modest numbers of antennas and relays, e.g. > 4 and 
K >2. Further improvements can be achieved by a process we 
term RS which addresses the possibility of mismatching poor 
first phase channels with optimized second phase channels as 
well as reducing the cardinality of fix- 

By removing one or more relays based on their MSB 
performance from consideration by (Il4l . fix can be optimized 



and its cardinality improved without overly restricting the 
second-phase channels available to the TDS process. The 
selection of the single highest MSB relay can be expressed 
as a discrete maximization problem given by 

jopt = arg max r^r^ , ] 



K 

arg max^^^E - Wj^k[i]rsr^\ 



(15) 



k=l 



where JIr is the set of candidate relays and Wj^k[i] is the 
MMSB filter for the k^^ symbol at the j^^ relay. On the 
solution of ( fTSl ), a refined subset, CIt G ^It, is generated by 
removing members of VIt which involve transmission from 
relay jopt, i.e. members of fiy where [a],, d - ■ ■'^r d\^^- 

3 opt 3 opt 

TDS then operates with this subset, where \CtT\ = (■^^'■~^')- 
Bxtension to the selection of multiple relays involves summing 
the MSB from candidate relays and populating VIb. with sets 
of these relays. However, the selection of the number of 
relays to remove is vital, as too high a value will result in 
a overly restricting the second phase channels available to the 
TDS process therefore increasing the probability of a channel 
mismatch. 



C. Cooperative MMSE Channel Estimation 

The next task that is necessary for the interference mit- 
igation and resource allocation is to compute the channel 
gains of the links of the cooperative system. In order to 
estimate the channel in the cooperative system under study, 
let us first consider the transmitted signal for user k, Xk[i\ = 
Bk[i]Ak[i]Ckhk[{\ = QjJ\i\hk\i], and the covariance matrix 
given by 

i?= [r[z]r^[z]] 

K 

k=l 
K 



k=l 

(16) 

A linear estimator of hk[i\ applied to r[i\ can be represented 
as hk[i] = Fj^r[i]. The linear MMSB channel estimation 
problem for the cooperative system under consideration is 
formulated as 

Fk.opt = argmin£;[||/ifc[i] - /ife[i]|P] 

(17) 

= argmin£;[||/ife[i] - Tf r[i]|p]. 

Tk 

Computing the gradient terms of the argument and equating 
them to zero yields the MMSB solution 

Fk,opt = R 'Pk, (18) 
E[r[i]h^[i]] = Q,[i]E[h,[i]h^[L]] = 



where P 



k 
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Qf,[i]Phf.. Using the relation hk[i] ~ Fj^r[i], we obtain 

K 

k=l 

(19) 

The expressions in (fT9] l require matrix inversions with cu- 
bic complexity ( 0{J'^)), however, this matrix inversion is 
common to (T3[ and needs to be computed only once for 
both expressions. In what follows, computationally efficient 
algorithms with quadratic complexity (0(J^)) based on an 
alternating optimization strategy will be detailed. 

IV. Adaptive Algorithms 

In this section, we present algorithms to compute the pa- 
rameters of interest and the expressions derived in the previous 
section with lower computational complexity. Specifically, we 
develop adaptive RALS algorithms using a method to build 
the group of G users based on the power levels, and then 
we employ an alternating optimization strategy for efficiently 
estimating the parameters of the receive filters, the power allo- 
cation vectors and the channels. Despite the joint optimization 
that is associated with a non-convex problem, the proposed 
RALS algorithms have been extensively tested and have not 
presented problems with local minima. 

A. Group Allocation and Channel Estimation 

The first step in the proposed strategy is to build the group of 
G users that will be used for the power allocation and receive 
filter design. A RAKE receiver Q, which is equivalent to a 
filter matched to the signature sequence of the desired user or 
the spatial signature of the desired data stream will be used 
for the group allocation. The RAKE receiver is employed to 
obtain zf^'^^[i\ = {C khk\i])" r\i] = p"[i\r\i]. The group is 
then formed according to 

compute the G largest jz^'^^WI, k = l,2,...,K. 

(20) 

The design of the RAKE and the other tasks require channel 
estimation. The power allocation, receive filter design and 
channel estimation expressions given in (fTSI i. ( fTsl l and 
respectively, are solved by replacing the expected values with 
time averages, and RLS-type algorithms with an alternating 
optimization strategy. In order to solve ([19) efficiently, we 
develop a variant of the RLS algorithm that is described by 



I r\i\ 



(21) 



where <5j,[i] = Bk[i]Ak[i]Ck, the estimate of the inverse of 
the covariance matrix R [i] is computed with the matrix 
inversion lemma ISTl 



'^R[i - l]r[i] 



(22) 



l + a~'^rfi[i\R[i - l]r[i\' 
R[i\ a-'^R[i - 1] - a~^k[i\r"[i\R[i - 1], (23) 



and 



where a is a forgetting factor that should be close to but less 
than 1. 

B. Joint Interference Suppression and Power Allocation 

The approach for allocating the power within a group is 
to drop the constraint, estimate the quantities of interest and 
then impose the constraint via a subsequent normalization. The 
group-based power allocation algorithm is computed by 

as,k[^ = Rs,k['i]Ps,k['^ 

= Rs.kmaps^kli - 1] + bk[i]vk[i]) (25) 
= as.k[i - 1] + ^a[i\ks,k[^\^ 

where S,a[i] ~ bk[i] — 6,5 ki''- ^ M'^kii] is the a priori error, 
Vk[i] = [i]Pg [i]wk[i] is the input signal to the recursion 

a^^Rs,k['i - l]^fcW 



k 



(26) 



1 + a-^v^[i\Rs,k[^ - '^]vk[i\ 

Rs.k[i] =- a~^Rs.k[^ - 1] - a-^ks.k[i]Vk[^Rs.k[^ " !]■ 

(27) 

The normalization a5_fc[i] <— Pq (is.k[i\/\\0'S,k[''']\\ 1^ then 
performed to ensure the power constraint. 
The linear receive filter is computed by 



Wk\i] = Wk[i - 1] + k[i\C\i], 



(28) 



where the a priori error is given by ^[i] = bk[i] — w^[i 
l]r[i] and k[i] is given by ( |22] |. The proposed scheme employs 
the algorithm in ( l20l i to allocate the users in the group and 
the channel estimation approach of (l2ni-(l24li. The alternating 
optimization strategy uses the recursions (|25t and (128) with 
1 or 2 iterations per symbol i. 

C. Transmit Diversity Selection and Relay Selection Based on 
Discrete Stochastic Gradient Algorithms 

In this part, we describe a low-complexity solution to the 
joint TDS and RS problem based on a discrete stochastic 
gradient algorithm (DSA) that can compute the optimal com- 
binatorial solution outline in (fT4l i and ( fTSl l with a substantially 
reduced cost as compared with the exhaustive search. We 
present a pair of low-complexity DSA that was first reported in 
1331 . II52I . which jointly optimizes RS and TDS in accordance 
with ( fT4b and ( fTSl ), and converges to the optimal exhaustive 
solution. 

The RS portion of the DSA is given by the algorithm of 
Table H] At each iteration the MSE of a randomly chosen 
candidate relay {j'^) (step 2) and that of the worst performing 
relay currently known (j^) are calculated (step 3). Via a 
comparison, the higher MSE relay is designated for the 
next iteration (step 3). The current solution and the relay 
chosen for removal {j) is denoted as the current optimum and 
is the relay which has occupied most frequently over the 
course of the packet up to the i*^ time instant; effectively 
an average of the occupiers of . This averaging/selection 
process is performed by allocating each member of fij^ a 
\^r\ X 1 unit vector, v;, which has a one in its corresponding 
position in 51 j^, i.e., v^h [i] is the label of the worst performing 
relay at the iteration. The current optimum is then chosen 
and tracked by means of a jilit'l x 1 state occupation probability 



Ph, [1] = aPn, [z - 1] + hk[i - 1]K - 1] 



(24) 
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TABLE I 

Proposed discrete stochastic joint TDS and RS algorithm 



1, 7rfl[l,j] = 



Step 

1. Initialization 

choose j[l] G UrJ^II] e Ur, 7TR[l,j[l 
for J 

2. For the time index i = 1,2, N 
choose 6 

3. Comparison and update of the worst performing relay 

ifT[i,r,r^Ci,-i] > [i= then jM'Ii + l] = 

otherwise j'^I* + 1] = 

4. State occupation probability (SOP) vector update 

■^R[i + 1] = TflH + mVjW[i+i] - TTRli]) where fi[i] = 1/i 

5. Determine largest SOP vector element and select the optimum relay 
if 7TR [i + + 1]] > 77R,[i + then j[i + 1] = + 1] 

otherwise j [i + 1] = j [i] 

6. TDS Set Reduction 

remove members of Qt which utiHze relay j[i + 1] (Qt ^ f^r) 



(SOP) vector, tvj^. This vector is updated at each iteration 
by adding Vjw [i + 1] and subtracting the previous value of 
TTfl (step 4). The current optimum is then determined by 
selecting the largest element in tt^ and its corresponding entry 
in flu (step 5). Through this process, the current optimum 
converges towards and tracks the exhaustive solution ll52l . An 
alternative interpretation of the proposed algorithm is to view 



the transitions, j 



+ as a Markov chain and the 



members of VLji as the possible transition states. The current 
optimum can then be defined as the most visited state. 

Once RS is complete at each time instant, set reduction 
(r^T — > step 6) and TDS can take place. To perform TDS, 
modified versions of steps 1 — 5 are used. The considered set 
is replaced, Q,r — s- VIt; the structure of interest is replaced, 
j — > Tr\ the best performing matrix is sought — > 7~f ; 
the SOP vector is replaced tt^j — > Ttf and C ^ T from 
(fT4l i. Finally, the inequality of step 3 is reversed to enable 
convergence to the lowest MSB TDS matrix which is then 
feedback to the relays in the form of 1-bit per relay antenna. 

Significant complexity savings result from the proposed 
algorithm; savings which increase with K, rir and the number 
relays removed in the RS process. For example, when rir = 10, 
K = 2 and 4 relays are removed, the number of complex 
multiplications for MMSE reception and exhaustive TDS, 
exhaustive TDS with RS, iterative TDS and iterative TDS 
with RS are 5.8 x 10**, 1.7 x 10^ 1.8 x 10^ and 5.9 x 10'^, 
respectively, for each time instant. 

V. Analysis and Requirements of the Algorithms 

In this section, we assess the requirements of the proposed 
and existing algorithms for cross-layer design in terms of 
computational complexity and number of feedback bits. The 
basic idea is to show the computational cost of the algorithms 
presented and compare them with those of existing techniques 
for interference mitigation and/or resource allocation. 

A. Computational Complexity Requirements 

We discuss here the computational complexity of the pro- 
posed and existing algorithms. Specifically, we will detail the 
required number of complex additions and multiplications of 
the proposed JPAIS-GBC algorithms and compare them with 



interference suppression schemes without cooperation (NCIS) 
and with cooperation (CIS) using an equal power allocation 
across the relays. Both uplink and downlink scenarios are 
considered in the analysis. In Table I we show the computa- 
tional complexity required by each recursion associated with 
a parameter vector/matrix for the JPAIS-GBC with G ~ K, 
which is more suitable for the uplink. 

table II 

Computational complexity of algorithms with a global 
power constraint g ^ k . 

Number of operations per symbol 



Parameter 


Additions 


Multiplications 


W\i] 


2[jf 




+2K{J) 


+2K{J) 




-J+l 


+3J + 1 




3K{K{nr + 1)) 


K{K{nr + l)f 




+K{nr + 1){L-1) 


+4:{K{nr + l)f 


OLT [i] 


+K{M{nr + 1)) 


+ {K + L){K{nr + l)f 
-{K{nr + l)f 




+K{K{nr + 1)) 
+6{K(nr + l)f 




+K{MQ) 




+ 3K{nr + 1) +1lr+2 


+nr 






+^(K{nr + If 




+5KQ 


+6KQ 




+3 


+1 



In Table II we show the computational complexity required 
by each recursion associated with a parameter vector for the 
JPAIS-GBC algorithm, which is suitable for both the uplink 
and the downlink. A noticeable difference between the JPAIS- 
GBC with G = K and G = 1 is that the latter is employed 
for each user, whereas the former is used for all the K users 
in the system. Since the computation of the inverse of R[i] is 
common to all users for the uplink in our system, the JPAIS- 
GBC with G = is more efficient than the JPAIS-GBC with 
G = 1 computed for all the K users. 

table iii 

Computational complexity of algorithms with 
individual power constraints (g = 1). 



Number of operations per symbol 



Parameter 


Additions 


Multiplications 




2{jr 




Wk[i] 


+j 


+5J 




+1 


+1 




2{nr + lf 


3(n, + 1)2 




+3{nr + 1) 


+7(n,. + 1) 


at [i] 


+ JL 


+JL 




+Q 


+Q 




-3 


+3 




mr 


GiQf 


hk[i\ 


+5MQ 


+MQ 




-5(n, + 1) + 3 


+4(?ir + 1) + 1 



The recursions employed for the proposed JPAIS-GBC with 
G = K and the JPAIS-GBC with G = 1 are general and 
parts of them are used in the existing CIS and NIS algorithms. 
Therefore, we can use them to describe the required compu- 
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tational complexity of the existing algorithms. In Table III we 
show the required recursions for the proposed and existing 
algorithms, whose complexity is detailed in Tables I and II. 

TABLE IV 

Computational complexity of the proposed JPAIS and 

EXISTING algorithms. 



Algorithm 


Recursions 






JPAIS-GPC(G = K) (Uplink) 


W[i\ 




M 


M 


JPAIS-GBC (G = 1) (Downlink) 


Wk [i] 


, ak [i 


,h 


k\i] 


CIS (Uplink) 


W[i\ 


a,T[i 


is 


fixed 


CIS (Downlink) 


Wk.[i] 


• ak[i 


is 


fixed 


NCIS (Uplink) 


W[i\ 


with 


rir 


= 


NCIS (Downlink) 


Wk[i] 


with 




= 



G = 1 . A typical number of bits required to quantize each 
coefficient of the vectors and via scalar quantization 
is n;, = 4 bits. More efficient schemes employing vector 
quantization ll53l . 1541 and that take into account correlations 
between the coefficients are also possible. 

For the uplink (or multiple-access channel), the base station 
(or access point) needs to feedback the power levels across the 
links to the K destination users in the system. With the JPAIS- 
GBC with G ~ K algorithm, the parameter vector with 
{rir + l)Knb bits/packet must be broadcasted to the K users. 
For the JPAIS-GBC algorithm with G = 1, a parameter vector 
ak with {rir + l)nb bits/packet must be broadcasted to each 
user in the systems. In terms of feedback, the JPAIS-GBC 
algorithm with G = 1 is more flexible and may require less 
feedback bits if there is no need for a constant update of the 
power levels for all K users. 
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Fig. 3. Computational complexity in tenns of the number of complex 
multiplications of the proposed and existing schemes for the uplink. 

In Fig. [3j we illustrate the required computational complex- 
ity for the proposed and existing schemes for different number 
of relays (?ir)- The curves show that the proposed JPAIS-GBC 
with G = K and JPAIS-GBC with G = 1 are more complex 
than the CIS scheme and the NCIS. This is due to the fact 
that the power allocation and channel estimation recursions 
are employed. However, we will show in the next section 
that this additional required complexity (which is modest) can 
significantly improve the performance of the system. 

B. Feedback Channel Requirements 

The JPAIS algorithms presented so far for cross-layer design 
require feedback signalling in order to allocate the power 
levels across the relays. In order to illustrate how these 
requirements are addressed, we can refer to Fig. H which 
depicts the structure for both the data and feedback packets. 
The data packet comprises a number of allocated bits for 
training (A^ti), for synchronization and control (A'sync) and 
the transmitted data (iVdata)- The feedback packet requires the 
transmission of the power allocation vector ax for the case of 
the JPAIS-GBC algorithm with G = K, whereas it requires 
the transmission of for each user for JPAIS-GBC with 



Data Packet Structure 



Feedback Packet Structure 
(JPAIS-GBC wtilr G=K) 



Feedback Packet Structure 
(.IPAIS-GBC with G=l) 



{n, + l)Knj, bits/feedback packet 



A" times (tIj. + l)Knb bits/feedback packet 



Fig. 4. Proposed structure of the data and feedback packets. 

For the downlink (or broadcast channel), the K users must 
feedback the power levels across the links to the base station. 
With the JPAIS-GBC algorithm with G = K, the parameter 
vector ut with (n^ + l)Knt bits/packet must be computed 
by each user and transmitted to the base station, which uses 
the aT vector coming from each user An algorithm for data 
fusion or a simple averaging procedure can be used. For the 
JPAIS-GBC algorithm with G = 1, a parameter vector 
with [n,. + l)ni, bits/packet must be transmitted from each 
user to the base station. In terms of feedback, the JPAIS-GBC 
algorithm with G = 1 requires significantly less feedback bits 
than the JPAIS-GBC with G = A' in this scenario. 

The MIMO TDS and RS scheme can be interpreted as a 
1-bit power allocation scheme and therefore achieves perfor- 
mance improvements whilst utilizing the minimum number 
of feedback bits per antenna. Consequently, the feedback 
requirements per update of a cooperative MIMO system using 
TDS and RS is given by the total number of relay transmit 
antennas UrK. This minimal feedback allows optimization of 
the system whilst maintaining the capacity of the system with 
regards to the transmission of useful data. 

VI. Simulations 

In this section, we illustrate with Monte-Carlo simulations 
the performance of the cross-layer algorithms described in this 
chapter Specifically, we assess the performance in terms of 
the bit error ratio (BER) of the JPAIS scheme and adaptive 
algorithms with group-based power constraints (GBC). The 
JPAIS scheme and algorithms are compared with schemes 
without cooperation (NCIS) and with cooperation (CIS) ll26l 
using an equal power allocation across the relays (the power 
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allocation in the JPAIS scheme is disabled). We also assess 
the proposed algorithms for transmit diversity selection and 
relay selection (Iterative TDS with RS) are presented and 
comparisons drawn against the optimal exhaustive solutions 
(Exhaustive TDS with RS), the unmodified system (No TDS), 
and the direct transmission (Non-Cooperative) . 

A. DS-CDMA System 

A DS-CDMA network with randomly generated spreading 
codes and a processing gain = 16 is considered. The 
fading channels are generated considering a random power 
delay profile with gains taken from a complex Gaussian 
variable with unit variance and mean zero, L = 5 paths 
spaced by one chip, and are normalized for unit power The 
power constraint parameter Pa fc is set for each user so that 
the designer can control the SNR (SNR = P^./t/f^) and 
Pt = Pg + {K — G)PA.k^ whereas it follows a log-normal 
distribution for the users with associated standard deviation 
equal to 3 dB. The DF cooperative protocol is adopted and 
all the relays and the destination terminal use either linear 
MMSE, which have full channel and noise variance knowl- 
edge, or adaptive receivers. The receivers are adjusted with 
the proposed RALS with 2 iterations for the JPAIS scheme, 
and with RLS algorithms for the NCIS and CIS schemes. 
We employ packets with 1500 QPSK symbols and average 
the curves over 1000 runs. For the adaptive receivers, we 
provide training sequences with A'ti = 200 symbols placed 
at the preamble of the packets. After the training sequence, 
the adaptive receivers are switched to decision-directed mode. 

The first experiment depicted in Fig. |5] shows the BER 
performance of the proposed JPAIS scheme and algorithms 
against the NCIS and CIS schemes with = 2 relays. The 
JPAIS scheme is considered with the group-based power con- 
straints (JPAIS-GBC). All techniques employ MMSE or RLS- 
type algorithms for estimation of the channels, the receive 
filters and the power allocation for each user The results show 
that as the group size G is increased the proposed JPAIS 
scheme and algorithms converge to approximately the same 
level of the cooperative JPAIS-MMSE scheme reported in ||28l , 
which employs G = K for power allocation, and has full 
knowledge of the channel and the noise variance. 

The proposed JPAIS-GBC scheme is then compared with a 
non-cooperative approach (NCIS) and a cooperative scheme 
with equal power allocation (CIS) across the relays for 
Ur = 1,2 relays. The results shown in Fig. |6] illustrate the 
performance improvement achieved by the JPAIS scheme and 
algorithms, which significantly outperform the CIS and the 
NCIS techniques. As the number of relays is increased so 
is the performance, reflecting the exploitation of the spatial 
diversity. In the scenario studied, the proposed JPAIS-GBC 
with = 3 can accommodate up to 3 more users as compared 
to the CIS scheme and double the capacity as compared with 
the NCIS for the same BER performance. The curves indicate 
that the GBC for power allocation with only a few users is 
able to attain a performance close to the JPAIS-GBC with 
G = K users, while requiring a lower complexity and less 
network signalling. A comprehensive study of the signalling 
requirements will be considered in a future work. 

The next experiment considers the average BER perfor- 
mance against the normalized fading rate /^T (cycles/symbol), 
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Fig. 5. BER performance versus number of symbols. Parameters: AF proto- 
col, At = Afc = 0.025 (for MMSE schemes), a = 0.998, Rs]k[i\ = O-d-'' 
and PC^[i] = 0.01/. 
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Fig. 6. BER performance versus SNR and number of users for the optimal 
linear MMSE detectors. Pai'ameters: AF protocol, a = 0.998, Rq = 
0.01/ and k~^[i\ = 0.01/. 



as depicted in Fig. [7] The idea is to illustrate a situation where 
the channel changes within a packet and the system transmits 
the power allocation vectors computed by the proposed JPAIS 
algorithms via a feedback channel. In this scenario, the JPAIS 
algorithms compute the parameters of the receiver and the 
power allocation vector, which is transmitted only once to the 
mobile users. This leads to a situation in which the power 
allocation becomes outdated. The results show that the gains of 
the proposed JPAIS algorithms decrease gradually as the fdT 
is increased to the BER level of the existing CIS algorithms for 
both n,. = 2 and ti^ = 4 relays, indicating that the power al- 
location is no longer able to provide performance advantages. 
This problem requires the deployment of a frequent update of 
the power allocation via feedback channels. 

The last experiment, shown in Fig. [8] illustrates the averaged 
BER performance versus the percentage of errors in the 
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Fig. 7. BER performance versus /^T for the AF protocol. The parameters 
of the adaptive algorithms are optimized for each /^T. 
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Fig. 8. BER performance versus percentage of feedback errors. 



feedback channel for an uplink scenario. Specifically, the 
feedback packet structure is employed and each coefficient 
is quantized with 4 bits. Each feedback packet is constructed 
with a sequence of binary data (Os and Is) and is transmitted 
over a binary symmetric channel (BSC) with an associated 
probability of error Pg- We then evaluate the BER of the 
proposed JPAIS and the existing algorithms against several 
values of the P^. The results show that the proposed JPAIS 
algorithms obtain significant gains over the existing CIS 
algorithm for values of < 0.1%. As we increase the rate 
of feedback errors, the performance of the proposed JPAIS 
becomes worse than the CIS algorithms. This suggests the use 
of error-control coding techniques to keep the level of errors 
in the feedback channel below a certain value. 



the unmodified system where all antennas are active (No TDS), 
and the direct transmission (Non-Cooperative). Plots of the 
schemes with TDS only (Exhaustive TDS, Iterative TDS) 
are also included to illustrate the performance improvement 
obtained by RS. Equal power allocation is maintained in each 
phase so that the total transmit bit power of the relays is 
unity. RLS channel estimation (CE) is used where all auxiliary 
matrices are initialized as identity matrices and estimation 
matrices are zero matrices, and the exponential forgetting 
factor is 0.9. Each simulation is averaged over 1000 packets 
(A^p), each with training sequences of 200 symbols. 
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Fig. 9. Cooperative DF MIMO BER performance versus the number of 
received symbols. 



Fig- 12] gives the BER convergence performance of the 
proposed algorithms. The iterative TDS with RS algorithm 
converges to the optimal BER as does TDS with RS and CE, 
albeit in a delayed fashion due to the CE. The TDS with RS 
scheme exhibits quicker convergence and lower steady state 
BER. These results and the interdependence between elements 
of the algorithm confirm that both the RS and TDS portions 
of the algorithm converge to their exhaustive solutions. 

Fig. [To] shows the BER versus SNR performance of the 
proposed and conventional algorithms. Increased diversity 
has been achieved without sacrificing multiplexing gain and 
illustrates that although the maximum available diversity ad- 
vantage decreases from Af (n^ + 1) to M{Ksub/K + 1) with 
RS with TDS because fewer antennas are active, the actual 
diversity achieved has increased. These diversity effects can 
be attributed to the removal of poor paths and therefore a 
lower probability of first phase and second phase channel 
mismatch but also the increase in transmit power over the 
remaining paths. The largest gains in diversity are present 
in the 15 — 25dB region and begin to diminish above this 
region because relay decoding becomes increasingly reliable 
and lower power paths become more viable for transmission. 



B. MIMO System 

In this part, simulations of the proposed algorithms (Iterative 
TDS with RS) are presented and comparisons drawn against 
the optimal exhaustive solutions (Exhaustive TDS with RS), 



VII. Extensions and Suggestions for Future Work 

The algorithms for joint resource allocation and interference 
mitigation described in this chapter are quite general and 
can be employed in a variety of wireless communication 
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Fig. 10. Cooperative FD MIMO BER perfomiance versus SNR. 



systems that are equipped with cooperative techniques. These 
include orthogonal-frequency-division-muhiplexing (OFDM) 
ifSTl . single-carrier systems with frequency-domain equalisa- 
tion (SC-FDE) 1581 and ultra- wide band (UWB) systems 1591 . 

Possible extensions include the incorporation of more ad- 
vanced interference mitigation strategies than linear schemes. 
These include nonlinear detection techniques such as succes- 
sive interference cancellation ifTTI . ll38l . decision feedback 
strategies l37l . Il39l and sphere decoders ||23l . The detection 
algorithms could also be considered with space-time coding 
schemes Q, IfTOl . channel coding and iterative processing 
approaches ||39l . l62l . 

Another complementary set of techniques comprises algo- 
rithms for adaptive parameter estimation. These methods are 
fundamental to estimate key parameters such as channel gains, 
amplitudes and receive filters, whilst keeping the complexity 
low and being able to track variations of the environment. 
Amongst the adaptive parameter estimation techniques, a 
designer can choose between supervised and blind approaches 
ISTl . Bhnd techniques l63l - l68l are appealing as they can 
increase the spectral efficiency of wireless systems. This is 
especially relevant for cooperative systems as they require 
extra signalling for cross-layer design. Supervised adaptive 
algorithms usually rely on training sequences that are sent at 
the beginning of each data packet ||69l . ifTOl . One fundamental 
issue in the choice of the adaptive parameter estimation 
algorithm is the speed of convergence and the tracking perfor- 
mance. The literature suggests that reduced-rank algorithms 
l43l - l49l . I72l-l74l are very atti-active choices when fast 
training and accurate tracking are important issues. 

VIII. Concluding Remarks 

We have presented in this work joint iterative power allo- 
cation and interference mitigation techniques for DS-CDMA 
and MIMO networks which employ multiple hops and the AF 
and DF cooperation protocols. A joint constrained optimiza- 
tion framework and algorithms that consider the allocation 
of power levels across the relays subject to group power 



constraints and the design of linear receivers for interfer- 
ence suppression were proposed. A scheme for joint transmit 
diversity optimisation and relay selection along with linear 
interference suppression has also been detailed and applied to 
MIMO systems. A study of the requirements of the proposed 
and existing algorithms in terms of computational complexity 
and feedback channels has also been conducted. The results of 
simulations have shown that the proposed algorithms obtain 
significant gains in performance and capacity over existing 
non-cooperative and cooperative schemes. 
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